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We propose the use of quantum optical systems to perform universal simulation of quantum 
dynamics. Two specific implementations that require present technology are put forward for il- 
lustrative purposes. The first scheme consists of neutral atoms stored in optical lattices, while 
the second scheme consists of ions stored in an array of micro-traps. Each atom (ion) supports 
a two-level system, on which local unitary operations can be performed through a laser beam. A 
raw interaction between neighboring two-level systems is achieved by conditionally displacing the 
corresponding atoms (ions). Then, average Hamiltonian techniques are used to achieve evolutions 
in time according to a large class of Hamiltonians. 



PACS numbers: 03.67.-a, 03.65.-w, 32.80.Pj, 42.50.-p 

I. INTRODUCTION 

Simulating quantum systems on a classical computer is 
known to be hard. Consider a set of two-level quantum 
systems, say N spin-1/2 particles placed at the sites of 
some regular lattice, that interact with each other. The 
number of parameters required to describe the state of 
these spins grows exponentially with N jl], |^|. For in- 
stance, the state of N = 50 spin-1/2 systems is specified 
by 2 50 w 10 15 numbers, whereas a 2 50 x 2 50 matrix, i.e. 
with ps 10 30 entries, needs to be exponentiated in order to 
compute its time evolution. Therefore, a device process- 
ing information according to the laws of classical physics 
— and in particular any present computer — is unable to 
efficiently simulate the dynamics of these spins when N 
becomes large. 

This fact has severe consequences in the study of 
condensed matter systems. Our understanding of a 
large spectrum of collective quantum phenomena, in- 
cluding ferromagnetism, antiferromagnetism, conductor- 
insulator transitions, superconductivity and quantum 
Hall effect, is dramatically hindered by the impossibil- 
ity of tracing quantum dynamics. Suppose the N spins, 
whose 2 JV -dimensional Hilbert space H decomposes as 

H N =H (1) ®H {2) ®---®H {N \ H {k) =C 2 , (1) 

display a property that we want to study, originating, 
say, in a given collective effect, and we conjecture that 
some Hamiltonian Hn acting on Hjv is able to account 
for this effect. Since a classical computer can not effi- 
ciently simulate a quantum evolution according to -ff/v, 
the question whether Hn successfully describes the prop- 
erty of interest can in general not be answered. And 
thus, for instance, the simple Hubbard model, believed 
to explain a wide range of electromagnetic properties of 



condensed matter systems -high temperature supercon- 
ductivity among them-, remains unsolved after decades 
of study §|. 

As a matter of fact, Feynman's initial motivation for 
constructing a quantum computer was the efficient sim- 
ulation of quantum dynamics 0] . Building on the obser- 
vation that a quantum system can be used to simulate 
another quantum system, he conjectured the existence of 
a universal quantum simulator UQS. Since then, several 
authors have analyzed this possibility ||. A UQS is a 
controlled device that, operating itself at the quantum 
level, efficiently reproduces the dynamics of any other 
many-particle quantum system that evolves according to 
short range interactions. Here, the assumption of some 
degree of locality in the interactions, implying that the 
multi-particle Hamiltonian Hjy = Hi is a sum of 
terms Hi each one involving only a few neighboring sys- 
tems, is important to achieve an efficient simulation. In 
most cases of interest this requirement happens to be ful- 
filled. Consequently, a UQS could be used to efficiently 
simulate the dynamics of a generic many-body quantum 
system and in this way function as a fundamental tool 
for research in quantum physics. 

On the other hand, the main present motivation for 
building a quantum computer comes from the expected 
exponential gain in efficiency of certain quantum algo- 
rithms with respect to their classical counterparts. Shor's 
efficient factorization of large numbers is so far the most 
celebrated milestone of quantum computation |6| . How- 
ever, for quantum computers to overcome classical ones 
in tasks such as factorization, they would have to coher- 
ently operate tens of thousands of two-level systems or 
quantum bits (qubits). This extraordinary enterprise re- 
quires technology that may only be at reach in several 
decades from now. Instead, simulating the dynamics of 
a few tens of qubits appears as a more feasible — and, 
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arguably, still immensely rewarding — task. 

In this paper we propose the use of quantum optical 
systems to realize a UQS. For illustrative purposes we put 
forward two specific schemes. First we consider the use 
of neutral atoms trapped in optical lattices R || . The 
second scheme consists of ions in an array of microtraps 
In both cases, some raw interaction is produced be- 
tween neighboring qubits by conditionally displacing the 
atoms (ions). Also, fast local unitary transformations are 
applied on the qubits, in a way that the evolution of the 
system effectively corresponds to that generated by some 
new, simulated Hamiltonian. In this way, a large class of 
Hamiltonians Hjq are produced. 

The present work aims, therefore, at bridging between 
previous theoretical developments concerning the simu- 
lation of quantum dynamics and its experimental realiza- 
tion. Among the main motivations for this contribution 
we encounter the belief that a universal quantum sim- 
ulator might be counted in the near future among the 
first demonstrated fundamental applications of quantum 
information science. 



II. OVERVIEW 

In this paper we propose, through two specific realiza- 
tions, the use of quantum optical systems to construct a 
UQS. That is, we put forward two specific schemes for 
the realization of a device that aims at simulating the 
dynamics of other quantum systems. 

We would like to note at the outset that from a com- 
putational point of view a UQS can be regarded as just 
a particular instance of a quantum computer and, con- 
versely, any quantum computer could be used to sim- 
ulate quantum systems. A simulation-oriented design, 
however, is likely to be decisive when it comes to actu- 
ally achieving quantum simulations with existing tech- 
nology |T(]] . It is this search for feasibility with presently 
available resources what possibly best defines the general 
philosophy of the present work. Instead of discussing fun- 
damental aspects of quantum dynamics simulation, here 
we will be primarily concerned with presenting both con- 
crete physical set-ups and protocols for its feasible imple- 
mentation. In particular, we will describe how to achieve 
the simulation of a considerably large class of interest- 
ing unitary dynamics that require particularly little time 
and control resources when performed in the proposed 
physical set-ups. 

A brief summary of the proposals, both having two- 
level systems or qubits as basic building blocks, is as 
follows. 

UQS1: Each qubit corresponds to two internal states 
|0) and |1) of a neutral atom. Two superposed, identical 
optical lattices are homogeneously filled with N atoms, 
one atom per lattice site. Local unitary transformations 
are performed on the qubits by shinning the atoms with a 
laser beam, whereas two-qubit interactions are produced 
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Figure 1: The simulation of a unitary evolution of N qubits 
(horizontal lines) according to a Hamiltonian Hn as in Eq. (^) 
is achieved by composing individual unitary transformations 
on each of the qubits (dark ellipses) with short gates that 
decompose into two-qubit components (empty ellipses). 



by displacing one of the optical lattices (which traps the 
atoms when their internal state is |1)) with respect to the 
other lattice (that traps the atoms when their internal 
state is |0)) [see figure |j. 

UQS2: Each qubit corresponds to two internal states 
|0) and |1) of an ion. N ions are stored in an array of 
microtraps, with one ion per trap. Again, local unitary 
transformations are performed on the qubits by means of 
a laser beam. In this scheme, an interaction between two 
selected qubits is achieved by pushing the corresponding 
ions with a force that only acts if the internal state of the 
ion is |1) [see figure [|. 

Schemes UQS1 and UQS2 share many features and we 
will take advantage of this fact by conducting a com- 
mon analysis in most of the paper. Both schemes are 
based on qubits, as in Eq. ([!]), described by the algebra 
of Pauli spin-1/2 operators. A second common feature 
concerns the way the multi-qubit system that constitutes 
the simulator can be externally addressed. The manipu- 
lation is divided into two classes of external interventions 
that produce, respectively, local unitary transformations 
on each qubit and two-body interactions between qubits. 
Simulations are achieved by combining these two possi- 
bilities (see figure Q). Through the interaction, short |ll|] 
two-qubit unitary transformations are produced. These 
are interspersed with fast one-qubit unitary transfor- 
mations, in a way that some average Hamiltonian ef- 
fectively guides the evolution of the qubits. Thus, we 
mimic decoupling and refocusing bang-bang techniques 
of nuclear magnetic resonance |l2| , although the present 
schemes also benefit notoriously from a more direct con- 
trol of the interactions. Of considerable use will be the 
characterization of non-local Hamiltonian simulation re- 
cently performed in the context of quantum information 

@ H- 

The basic features of these proposals, namely the use of 
spin-1/2 systems and the enforcement of one-qubit and 
two-qubit Hamiltonian evolutions, have practical impli- 
cations that are worth discussing in this introductory 
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part. By definition, a UQS is a device able to simulate ef- 
ficiently any multi-particle quantum system with range- 
restricted interactions and, as mentioned above, could in 
principle function as a general purpose quantum com- 
puter. However, technological limitations significantly 
reduce the tasks that it can accomplish in practice. This 
observation affects also the kind of systems that a specific 
implementation of a UQS can effectively simulate given 
some experimental possibilities. In particular, the very 
nature of our proposals makes them best suited for sim- 
ulating evolutions of systems whose building blocks are 
also two-level systems, and having a Hamiltonian 



Finally, in section VI we discuss how to improve pro- 



H N =Y / H {a) 



H (ab) 

a^b 



(2) 



that decomposes into one-qubit terms and two- 

qubit terms H^ ab \ [Here a and b are indices that label 
the qubits]. One finds, for instance, that for this class 
of systems the time required to perform a simulation us- 
ing our schemes is proportional to the time of simulated 
evolution, a very convenient fact Jl6| . Accordingly, here 
we will focus the discussion of simulation protocols pri- 
marily to the case of Hamiltonians of the form (||) . These 
include a variety of popular spin models, such as Ising in- 
teractions and both isotropic and anisotropic Heisenberg 
interactions j|. 

Finally, a brief comparison of the two proposals com- 
prises the following observations. What makes scheme 
UQS1 most attractive is its feasibility with present tech- 
nology, as supported by the results of current experi- 
ments with neutral atoms in optical lattices f| . Another 
appealing feature is the high degree of parallelism in the 
manipulation of the qubits. Its major drawback is that 
in present experiments only homogeneous local opera- 
tions can be performed — so that all qubits experience 
the same local unitary transformations — and two-qubit 
interactions necessarily follow a translationally invariant 
pattern. The degree to which each qubit can be indi- 
vidually addressed is an important issue, since it deter- 
mines whether the scheme constitutes a universal simu- 
lator. We will eventually suggest methods to operate on 
each qubit independently within scheme UQS1, but for 
most of the paper only homogeneous manipulations are 
assumed, implying that only a subclass of Hamiltonians 
(Q) can be simulated. Scheme UQS2, in turn, is distin- 
guished precisely by its higher degree of addressability. 
In this proposal each qubit can be individually manip- 
ulated, so that arbitrary (inhomogeneous) local unitary 
transformations can be performed. Also, the two-body 
interaction can be produced between arbitrary — but not 
too distant — pre-selected qubits, allowing for the simu- 
lation of any Hamiltonian of the form (||) . 

The rest of the paper is distributed as follows. Section 
III introduces some background material on the simula- 
tion of quantum dynamics. Section IV contains a dis- 
cussion of the schemes UQS1 and UQS2. In section [v] a 
series of examples are presented to illustrate how to sim- 
ulate the dynamics of N qubit systems in our set-ups. 



posal UQS1, and we consider, in the context of studying 
the ground state of a multi-qubit Hamiltonian, the effect 
of errors and imperfections in the simulation. 



III. SIMULATION OF QUANTUM DYNAMICS 

This section presents background material concerning 
the simulation of unitary quantum dynamics. We first 
describe how an effective Hamiltonian evolution can be 
obtained by time averaging over some other Hamiltonian 
evolutions. After a brief analysis of the time and control 
resources required in a simulation of this type, we re- 
view useful results of what is known among the quantum 
information community as non-local Hamiltonian sim- 
ulation using local unitary transformations [|l3| [l4[ [l5| , 
a technique that uses fast local unitary operations on 
qubits to effectively modify an existing interaction be- 
tween them. The reason for discussing this technique is 
that its requirements match the control possibilities of 
our two proposals. 

A starting observation concerning the simulation of 
quantum dynamics is that if a Hamiltonian K = X)i=i Ki 
decomposes into terms Ki acting on a small constant sub- 
space, then by Trotter formula jl7|, 



-iKr 



lim 

m — >oo 



-iKiT/rrig— iK? r /m 



-iK s r jra 



,(3) 

we can approximate an evolution according to Hamilto- 
nian K by a series of short evolutions according to each 
of the pieces Ki . Therefore we can simulate the evolution 
of the iV-qubit system of Eq. (|l|) according to Hamilto- 
nian Hn of Eq. (Q) by composing short one-qubit and 
two-qubit evolutions generated, respectively, by the one- 
qubit and two-qubit Hamiltonians and H^K 

In our schemes an evolution according to one-qubit 
Hamiltonians will be obtained directly by properly 
shinning a laser beam on the atoms or ions that host the 
qubits. Instead, two-qubit Hamiltonians H^ ah ^ will be 
achieved by processing some given interaction Hq that 
is externally enforced, as we explain in the following. 

Let us consider only two of the N qubits, that we de- 
note a and b. By alternating evolutions according to 
some available, switchable two-qubit interaction 
for times {ti} with local unitary transformations, one can 
achieve an evolution 



u(t) = n^ ex p(-^o (a6) *i)^ 

i=l 
n 

= Hexpf-iViH^vS), (4) 



where t = J2*i=o ^» = u ^ ® u i j with u i an d v i be- 
ing one-qubit unitary transformations, and the products 
are time-ordered. We proceed by considering the case 
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where the total time t is very small as compared to the 
time-scale of H^, so that — 1 <C tH^ <C 1. For nota- 
tional convenience, Hamiltonian H^ b ^ and time t are set 
to be dimensionless, with the eigenvalues of Hq being 
of order 1. Then, for t <C 1 we can approximate the ex- 
ponentials in Eq. by the first terms in their expansion 
in powers of t, to obtain 



U(t) ~ 1 



0(t 2 ), 



(5) 



where pi = U/t. That is, the evolution of the two qubits 
corresponds, effectively, to having the Hamiltonian 



H 



(ah) 
eff 



'mVH^V? +(D{t), 



(6) 



1=1 



acting on them for time t. Then, by concatenating several 
short gates U(t), 



U(t) 



0(t 2 ), 



(7) 



we can simulate -H^// f° r l ar g er times. 

Two parameters, the time cost c and the control com- 
plexity Xi can be used to characterize the above simu- 
lation. We consider the fast control limit, where local 
operations Vi are performed very fast as compared to the 
time scale of , so that the simulation time is deter- 
mined only by the use of H 1 ^^ . Suppose the aim is to 
simulate an evolution according to ij( ab ) = cH^f^ , c > 0, 
for a simulated time T". It follows from the linearity of 
Eq. (^) that this requires using for a time cT'. We 

define the time cost c as 



T 

rpi 1 



(8) 



so that it measures the time overhead required to simu- 
late H ( afc ) by iJ <j afc) . We remark that this measure of time 
resources implicitly assumes that the simulation time T is 
proportional to the simulated time T' . On the other 
hand, a simulation that takes time T is achieved by com- 
posing a number L of gates U(t) such that T = Lt. At 
each gate an error of order t 2 is introduced (see Eq. ([?[)), 
so that after time T the total error is of order e = Lt 2 , 
provided e C 1. For a given error e and simulated 
time T' , we find that the number of gates U(t) must 
be L = c 2 T' 2 /e. That is, in order to obtain a constant 
error e in a simulation for simulated time T', the total 
number L of gates U(t) must grow quadratically in T' , 
while the small time step t defining U(t) must decrease 
as t = e/(cT'). Then, if it takes n control operations Vi 
to perform each gate U(t), the total number of control 
operations — per unit of simulation time T — required to 
simulate H ( ab ) by Hq 1 is given by 



X 



nL 



ncT' 



(9) 



where \ defines the control complexity (per time unit) of 
the simulation. 

Summarizing, a two-qubit Hamiltonian and fast local 
operations allow us to simulate other two-qubit Hamil- 
tonians. The simulation time T is proportional to the 
simulated time T", and the required number \ of control 
operations per unit time grows linearly with T' ]l8( ] . The 
above analysis can now be carried out for N qubits with 
one-qubit and two-qubit Hamiltonians. A relevant as- 
pect is how the time cost c and the control complexity 
scale with N, for this may ultimately determine whether 
a given simulation is feasible. We find, for instance, that 
whenever parallel manipulation of all N qubits is possi- 
ble, the time cost c will not depend on N, whereas the 
control complexity will grow linearly in N, \ — ncNT' /e. 

We move now to consider the degree of independence 
between tt, and Vi in the control operation Vi = u\ ®u- 
of Eqs. (@)-@, which relates to the level of single-qubit 
addressability available in a given physical setting. The 
two extreme cases correspond to homogeneous manipula- 
tion, where the laser beam affects equally the two qubits, 
and thus u = v; and to inhomogeneous manipulation, 
where the laser beam can be sufficiently focused as to 
discriminate between the qubits and enforce independent 
evolutions on them, so that u and v can be arbitrary. In 
what follows we review the characterization, for homo- 
geneous and inhomogeneous manipulation, of the two- 
qubit Hamiltonian evolutions that can be simulated by 



Tj(ab) 



(10) 



for a real 7, — equivalently, by short gates 

- and by fast local operations. We quote results of 
]l4] , [ll|, which are time optimal, that is with minimal 
time cost c, and that happen to have a low — possibly 
optimal — control complexity \- As m Eq. (|l0|), we will 
express Hamiltonians in terms of the Pauli matrices 



1 

1 



-i 

1 



1 

-1 



(12) 



(i) Homogeneous manipulation |Tj]. By using homoge- 
neous local unitary transformations HLU as control op- 



erations V in Eq. 



that is V 



,.(«) 



<g> u 



(*>) 



simulate interaction Hamiltonians of the form 



H = 



E 



MijOi 



Mij = Mji, 



one can 



(13) 



which are symmetric under exchange of the qubits. As 
follows from the general analysis of Ref. UM , some ad- 
ditional restrictions on the real symmetric, 3x3 matrix 
M apply. 

Result 1. Hamiltonian H can be simulated by short 
gates according to in Eq. ( |l0| ) and fast HLU if and 
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only if the sign of 7 coincides with the sign of all non- 
vanishing eigenvalues of M. The time cost c of the 
simulation is c = (X^A*i)/7- 

Example 1. Consider the ferromagnetic [antiferromag- 
netic] Heisenberg interaction 



H = J 



(14) 



where J < [J > 0]. An evolution e~ lHT according to 
this Hamiltonian can be simulated by short gates (11), 
where 7 must be such that J7 > 0, alternated with fast 
HLU, according to the collection of n — 3 coefficients 
{pi} and local unitary {Vi} (see Eq, 



Pi 

P2 
P3 



i Vi = 1 ® 1, 



3' 



u 



1 — ia x 1 — ifia; 



3, v 3 = 



V2 
1 — icr 



V2 ' 



1 — ia 



y 



V2 



The time cost for simulating ( ^4[ ) with (|10|) is c = 3 J/7, 
and the control complexity x of Eq. (0) — or the number 
of control operations per unit time required in the simu- 
lation to obtain a small error e — is ncT'/e = 9JT' / (je). 

(ii) Inhomogeneous manipulation ]l5[ ]. The possibility 
to perform independent operations on each of the qubits 
translates into the ability to simulate all possible bipar- 
tite Hamiltonians. 

Result 2. Any two-qubit interaction Hamiltonian 



H = 



Mij0~i ® o~j 



(16) 



can be achieved by alternating the short two-qubit gate 
U of Eq. (|llj) with fast control operations of the form 
V = u(&v, that is inhomogeneous local unitary transfor- 
mations LU 15 1 . The time cost in time-optimal simula- 



tions is a function of the singular values {/ij > 0} of the 
real, 3x3 matrix M, and reads c = aWM- 

Example 2. An evolution according to the antisym- 
metric Hamiltonian 



H = J(a z (g) o y — G y <£> a z ), 



(17) 



can be achieved by composing gates U and a binary se- 
quence of LU characterized by 



1 1 + io~ T 

1 , 1 — io~ T 
P, - -, V 9 = _^-®1. 



(18) 



The time cost of the simulation is c = 2|J|/|7|, and the 
control complexity can be obtained from c and n = 2. 

We note here that inhomogeneous manipulation does 
not only have the advantage, as compared to homoge- 
neous manipulation, that it enables the simulation of 



arbitrary two-qubit Hamiltonians, but also that when 
dealing with N qubit systems, simulated interactions can 
depend on each couple of qubits. That is, unlike in the 
homogeneous case, the simulated Hamiltonian TJjv need 
not be translational invariant. A final comment concerns 
one qubit Hamiltonians of Eq. (||). In the homoge- 
neous case we find that must be the same for each 
qubit, whereas in the inhomogeneous case can be chosen 
independently for each qubit. 

With these basic techniques of Hamiltonian simulation 
in mind, we can now proceed to the remaining sections of 
the paper, that explain how to simulate iV-qubit quan- 
tum dynamics by using quantum optical systems. 



IV. PHYSICAL SET-UPS 

In this section we propose and discuss two physical set- 
ups that can be used to simulate many-body quantum 



(15) dynamics. In subsection IV A proposal UQS1, consisting 



of neutral atoms in optical lattices, is presented. In sub- 
we describe a second scheme, UQS2, that 



IV B 



section 

consists of ions trapped in an array of micro-traps. 



A. UQS1: Neutral atoms in optical lattices 

The first scheme we describe consists of neutral atoms 
trapped in optical lattices. Two internal states of the 
atoms, |0) and |1), define the qubits that are the relevant 
degrees of freedom for the simulated dynamics. We con- 
sider two identical, one-dimensional optical lattices with 
period d, each formed by the standing wave of two inter- 
fering laser beams that are superimposed in space p9[ . 
The lattices are such that they trap the neutral atoms de- 
pending on their internal state, |0) and |1). We assume 
that the double lattice is filled homogeneously with one 
atom per lattice site, and that the atoms are cooled to 
the vibrational ground state (this can be accomplished 
for instance by loading the lattice with a BEC (7[ [|). If 
required, the qubits can be initialized to state |0) through 
optical pumping. Their read-out can be achieved by 
measuring the phosphorescence of atoms after shining 
them with a proper laser. Depending on the degree of ad- 
dressability in the scheme, only collective measurements 
of the qubits may be feasible. 

Several proposals to perform coherent evolution of the 
joint state of two atoms trapped in optical lattices have 
been put forward |t], |2(J. In |21 , it was shown how to 
use optical lattice systems to simulate ferromagnetism 
and spin squeezing. We will base our scheme on the pro- 
posals and p^j , although some new techniques will be 
introduced. Note that one may also use different tech- 
niques to achieve a controlled interaction between neigh- 
boring atoms, including recent proposals to increase pos- 
sible gate fidelities of resulting phase gates by orders of 
magnitude p2]|. 
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Figure 2: Atoms in a double optical lattice. An interaction 
between adjacent qubits is achieved by displacing one of the 
lattices (which traps the atoms when their internal state is 
jl), light balls) with respect to the other (that traps the atoms 
when their internal state is |0), dark balls). In this way the |1) 
component of atom a approaches in space the |0) component 
of atom a + 1, and these collide in a controlled way. Then 
the two components of each atom are brought back together, 
and local unitary transformations can be enforced by shining 
a laser. By a sufficiently large, relative displacement of the 
two lattices, also interactions between more distant qubits 
can be achieved. Notice that the manipulation of the qubits 
is homogeneous and highly parallel. 



In the following we will explain how the interaction be- 
tween the internal degrees of freedom of different atoms 
can be accomplished. By adjusting the phases of the in- 
terfering laser beams, the two lattices are easily displaced 
with respect to each other, therefore displacing the atoms 
conditionally to their internal state. Performing this pro- 
cess adiabatically — such that the trapped atoms remain 
in the motional ground state — and choosing a relative 
displacement of one lattice period, we have that the |1) 
component of atom a meets the |0) component of atom 
a + 1, and the two atoms interact through controlled col- 
lisions for time t±. Then the lattices are returned (adia- 
batically) to their initial positions. In the absence of in- 
teractions, adiabaticity requires that \x(t)\ -C v sa where 
x is the relative position coordinate for the system of the 
two lattices, v osc w aou is the rms velocity of the atoms 
in the vibrational ground state, ui is the excitation fre- 
quency and ao is the size of the ground state of the trap 
potential 0. 

The result of the interaction process is that the 
|l)o|0)o+i ( as before, here a labels the position of the 
atoms) component of the wave functions picks up a cer- 
tain phase shift, which depends on the interaction time 
t\. This is a nonlocal evolution, generated by a Hamilto- 
nian whose interaction part reads 

^E^^g^, (19) 
a 

In addition, there are other phases originating in local 
Hamiltonians of the form ^ Q af^ ■ They can be removed 
by applying homogeneous local unitary operations on the 
atoms and will not be considered in what follows. There- 



fore, the resulting gate is of the form 

U x = e- l0lKl = Y[exp[-iOnr<?> (8 (20) 

a 

where 9\ can be adjusted by changing the interaction 
time t\. This gate, together with homogeneous local op- 
erations on the qubits, can be used to simulate evolutions 
according to other first neighbor interactions. 

Similarly, by displacing the lattices further, e.g. for 
j periods with respect to each other, and allowing for 
interaction time tj before returning to the initial position, 
each of the atoms interacts with its j th neighbor in a 
controlled way. In this case a gate of the form 

Uj = e~ ie i K i = Y[ expt-ifljO-W <g> eri Q+j) ], (21) 

a 

can be produced, where the Hamiltonian Kj is given by 
i^ = $>W®<7^), (22) 

a 

and 8j depends on the interaction time tj. In this way, 
two-qubit interactions between j th neighbors can also be 
simulated. 

Combining the above processes for different j by vis- 
iting all corresponding positions one can produce a total 
interaction Hamiltonian Hq which includes two-qubit in- 
teractions to first neighbors, second neighbors, etc. Note 
that part of the time required to produce these evolu- 
tions is spent in adiabatically shifting the lattice back 
and forth. For 9± = tt in Eq. (|2(]), the total time T re- 
quired for a complete cycle — taking adiabaticity require- 
ments into account — is at the order of few microseconds 
p3[ . This may be compared to a spontaneous emission 
lifetime of the single atoms in the lattice which ranges 
from several seconds to many minutes. 

This setup can be easily adapted also to 2D / 3D arrays 
of optical lattices. For example, nearest neighbor inter- 
action in a 2D square lattice is done by repeating the 
process for ID in both directions of the square lattice. 
Note that the whole process moves all the atoms at the 
same time, so that the interactions that are produced 
are invariant under translations. Homogeneous local op- 
erations on the qubits, the most feasible class of control 
operations on this setting, will not break this transla- 
tional symmetry. Thus, this scheme is mainly adequate 
to simulate translationally invariant Hamiltonians. 



B. UQS2: Ions in an array of micro-traps 

The basic architecture of this proposal consists of N 
ions, each one confined in an individual micro-trap ||, 
that we regard as a three-dimensional harmonic oscilla- 
tor. Thus we consider N micro-traps or independent 
three-dimensional harmonic oscillators, possibly gener- 
ated by electric or magnetic fields p4j and distributed 
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potential for ions A and -B reads 



Figure 3: Ions in an array of microtraps. An interaction 
between two qubits is achieved by conditionally pushing the 
corresponding ions with a force that only acts when their 
internal state is |1). The non-local evolution originates from 
a difference in the electrostatic potential affecting each ion, 
which depends simultaneously on the internal state of both 
ions. Local unitary operations are enforced on selected qubits 
by shinning a laser beam on the corresponding ions. 



in space according to some one-dimensional or two- 
dimensional pattern. Most of the present discussion will 
not depend on the particular geometry of the array. We 
only require the separation between neighboring traps to 
be such that a laser beam can address ions individually 
and that the Coulomb interaction between ions is not 
able to excite their vibrational state. 

The relevant degrees of freedom of the simulated dy- 
namics are two internal states of the ions, labeled as |0) 
and |1), that define a qubit. The initialization and read- 
out of these qubits has been discussed in [^5| . Since the 
separation between adjacent microtraps allows for indi- 
vidual addressing of the ions, local unitary operations 
on a given qubit can be performed by shinning a laser 
beam onto the corresponding ion, with laser frequency 
and phase chosen as not to affect its motion. 

Next we discuss how a two-qubit interaction can be 
achieved Suppose, for the time being, that the ions 
are in the ground state of the trap 1 and consider a force 
F(t) = y{t)huj/a \l)(l\ that acts on a selected ion only 
when its internal state is |1), producing a slight displace- 
ment in a particular direction. Here f(t) € [0, 1] is a 
smooth function of time, w is the frequency of the traps, 
ao = \ffi/2muo is the typical size of the ground state of 
the trap and m is the mass of the ion. It is important 
to note that this force can be chosen to act only on a 
particular set of ions in the array of micro-traps and that 
we can choose the direction of the displacement suffered 
by them. 

Suppose now the force F(t) acts on two ions, A and 
B, w hich suffer a conditional displacement in the direc- 
tion AE>. If we ignore the rest of ions in the array, the 



VAB = X) ^-[(ia-^)|l) Q (l|) 2 -5(<) 2 |l) a (l|] 

a=A,B 



(23) 



47re \d + x B - xaY 



where d is the distance between the centers of the traps, 
x(t) = f(t)ao and x a is the position operator for the ion 
a = A,B with respect to the original equilibrium posi- 
tion. Let Xa determine the minimum of this potential 
when the force is not acting, and let us replace the po- 
sition operators x a with the displacements around x^ 
and also redefine d — > d — Xj^ + x^ . We consider the 
case where \x a \ <C d, that allows for a perturbative ex- 
pansion of Coulomb potential; c\xaxb\/o^ <C 1, meaning 
that the interaction is a perturbation with respect to the 
traps: and e 2 -C 1, where e = . e , — \yz is the ratio 
of the Coulomb energy and the energy of one ion with 
respect to its neighboring trap. Under these conditions, 
if f(t) goes from to approximately 1 and then back to 
adiabatically (which requires \f(t)\ <C ui), the only effect 
on the ions will be the appearance of a phase 4>ij depend- 
ing on the internal state of the ions, i, j = 0, 1. 
The non-local content of this evolution is given by <j), 



m - 

„2 



dt 



Jo 



1 



1 



x B - xa 



Xb 



-nijj'e j dt XA{t)x B (t), 
'0 



d — xa d 
(24) 



where the contribution cf>ij has been computed by con- 
sidering the electrostatic potential at the time-dependent 
minimum (x A(t) , x b (t)) of Vab when the internal state 
of the ions is |i)^|j) B . 

The total effect of pushing the two ions back and forth 
is therefore to transform their internal states according 
to a unitary transformation U that, up to local terms 
that can be subsequently undone by applying a proper 
laser beam, reads 



U = e 



-iv cr Cw, 



e 2 a 2 1 ' T 



dtf{t) A f{t)i 



87TC0 d 3 Jo 

that is, according to the interaction Hamiltonian 



H, 



(AB) 



(25) 



(26) 



This condition will be relaxed later on. 



Then, arbitrary two-qubit Hamiltonian evolutions can be 
simulated between qubits A and B by alternating gates 
of the form ( p5|) and local operations on each qubit, as 
explained in the preceding section. 

Note that the conditional phases (fiij in ([m]) depend 
only on the mean position of the ions and not on the 
width of their wave functions. Hence they are insen- 
sitive to the temperature, and the ions need not be in 
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the ground state of the traps as the force starts acting. 
Notice also that both the external force F a (t) and the 
Coulomb potential due to the rest of ions may also con- 
tribute to the phases but they will not affect the 
non-local phase <f>. 

For later reference we point out that a force homo- 
geneously applied to all the ions naturally produces, be- 
cause of the 1/d 3 decay of 6 in Eq. ( |25| ) , a gate according 
to the Hamiltonian 



ff = Erf 1 3^ O) ®^ 



(6) 



a/6 



(27) 



where d a b denotes the distance between ions a and b. 

A final observation is that this natural l/|a — &| 3 de- 
cay in the strength of the interactions can be used to 
perform parallel processing on TV qubits, for large N. 
For instance, if we apply simultaneously the force F(t) 
to, say, qubits a and a + 1 and qubits a + 10 and a+11, 
the strength of the interaction between, say, qubit a and 
qubit a + 10 will be 10 -3 times that between qubits a 
and a+1, so that it may be neglected. In particular, this 
means that, for sufficiently large N, the time cost c of 
simulating a Hamiltonian Hn involving N qubits may 
be made not to depend on N. 



V. QUANTUM SIMULATOR AT WORK 

This section presents simulations that could be carried 
out with the proposals UQS1 and UQS2 of section [Tv| . 
Only a reduced number of simple examples have been 
selected for this discussion, that is not meant to be ex- 
haustive, but just intends to illustrate the broad range of 
applications that a universal quantum simulator would 
have in the study of condensed matter systems. 

We recall that in both physical set-ups of section IV 
one-qubit evolutions are achieved by shinning a laser on 
the systems (either atoms or ions) carrying the qubits, 
while two-qubit interactions are produced by condition- 
ally displacing these systems. Then, one can use the 
techniques of Hamiltonian simulation of section III to 
achieve evolutions of N qubits according to some effec- 
tive Hamiltonian H^. 

The first three examples we discuss correspond to 
Hamiltonians Hm that are invariant under translations. 
Their simulation does not require single-qubit address- 
ability, and therefore can be performed in any of the two 
set-ups UQS1 and UQS2. The fourth example requires 
the ability, inherent in proposal UQS2, of manipulating 
each qubit independently. 



and with dipole-dipole interactions. It is described by 
a Hamiltonian of the form 



1 



Hn = - 



J 



2 ~ld ab l 

a^tb 



» 



» 



(28) 



where er_ = |0)(1|, cr + = |1)(0|, and d a b denotes the 
distance between spins at lattice sites a and b. 

A possible interest to simulate Hamiltonian Hp comes 
from the fact that, in the 2D case, with anisotropic 
dipole interaction -that is, with the interaction strength 
J depending on the direction between spins- this simple 
model is known to have a spin-glass phase |2(|, which 
depends on the degree of anisotropy. 

An evolution according to Hp can be easily simulated 
if short evolutions according to a Hamiltonian 



Hn = 



a^b 



^^,d ab 3 



r («) 



(29) 



are enforced on the N qubits. In scheme UQS1, this can 
be accomplished by concatenating gates Uj of Eq. ( |2l| ) 
for different j, and by tuning the interaction times tj so 
that coefficients 8j decay with the cube of the distance 
between atoms. In scheme UQS2, Hamiltonian Hq can 
be produced by simply pushing all ions simultaneously 
with the conditional force F(t), as explained above (see 
Eq. (|27|)), since the 1/d 3 decay appears naturally there. 
In either case, short gates according to Hq can be alter- 
nated with homogeneous local unitary operations on the 
N qubits, as we did for N = 2 in section [II (see example 
1). If we consider the sequence {pi,Vi} of weights and 
local control operations given by 



1 T/ ft-ia x 

1 / 1 - i(7y 

P2 = r v >={— 



then, we have that the average Hamiltonian reads 



^pMHoV? 

i=l 



(30) 
(31) 



(32) 



a=jtb 



which is equivalent to Hp, as follows from the identities 
a x = {<7 + + cr_), (T y = -i(o+ - cr_). 



B. Ising and Heisenberg Hamiltonians 



A. Dipole-dipole interaction Hamiltonian 

We first consider a system of iV spin-1/2 particles 
placed at the sites of some ID or 2D regular lattice 



We move now to discuss the simulation of the Ising 
Hamiltonian 



Z ! 



(33) 



(a,b) 
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as well as that of the Heisenberg Hamiltonian 



H 



H 



E 

(a,b) 



r (o) 



r<«0 



r (a) 



(34) 

where the symbol (, ) means that the sums include only 
first neighboring sites in the array. 

Simulating these Hamiltonians appears as a promis- 
ing enterprise. For instance, the Ising model in a 2D 
triangular lattice and with J < can be used to study 
frustration effects, whereas the Heisenberg model in a 
2D triangular lattice has been discussed in the context 
of high temperature superconductivity and quantum Hall 
effect^. 

The simulation of Hi is particularly easy using scheme 
UQS1. In a ID setting, a short gate according to Hi 
can be achieved by displacing the lattices as explained 
in the previous section (recall Eq. (|l9|)). In 2D and 
3D, a lattice displacement in each spatial direction must 
be enforced. The achievement of only first-neighbor in- 
teractions in scheme UQS2 requires creating short gates 
according to the Hamiltonian where a and 

b are first-neighbor qubits. This is achieved by pushing 
the corresponding two ions. Each couple of ions has to 
be pushed at a different time, in order to avoid undesired 
interactions between ions that are not first neighbors. 
However, we already mentioned that the 1/d 3 decay in 
Eq. ( |25| ) allows for parallel processing of sufficiently dis- 
tant couples of ions. 

The Heisenberg model Hjj can be simulated from the 
Ising model Hi (see also pH). We have already described 
the N = 2 case in example 1 of section III. Similarly, in 



the general N case we can compose short gates according 
to Hi with homogeneous local operations according to 
the set {pi, Vi} given by 



Note that introducing local terms to a given simulated 
Hamiltonian can always be done, irrespectively of the 
Hamiltonian. Whether these local terms can depend on 
each qubit or must be homogeneous will depend on the 
available degree of individual addressability in the phys- 
ical scheme for simulation. 



Ising model with random coefficients and 
magnetic field 



Let us consider a Hamiltonian of the form 

H R = - \ Jab 4^ ® ^ + E B " a * ] > ( 38 ) 
(a,b) a 

that is, with first-neighbor Ising interaction, where coef- 
ficients J a b depend on the pair of qubits under consider- 
ation, and a qubit-dependent magnetic field along the x 
direction. 

In the ID case and with coefficients J a b randomly 
drawn from some probability distribution, this model has 
been used in the study spin glasses 26 and percolation 
It has also been considered in the context of solving 
classical problems with quantum annealing, as compared 
to classical annealing methods J2i| . 

Simulating this Hamiltonian requires single qubit ad- 
dressability. In the proposal UQS2, Hu can be simu- 
lated by producing short two-qubit gates U a b according 
to af^ eg) cri , where qubits a and b are first neighbors. 
One possibility is to assign a frequency v a b-, proportional 
to Jat, to each couple of qubits (a, b), and to perform 
gate U a b with that frequency. The magnetic field can 
be simulated by means of local operations, as we have 
described above. 



Pi 

P2 
P3 



3' 



Vi = l 



I v 2= ( 

3' 2 \ V2 

I v 3= ( IzJZl 

3' 3 { V2 



(35) 



In these two Hamiltonians one could add an extra term 
of a magnetic field 



(a) 



(36) 



where ft denotes a particular direction in real space. To 
simulate the presence of a magnetic field one has to apply 
a short unitary operation 



U = cxp [—i (a ■ n) St 



(37) 



where a = {c x ,(Ty,a z ), at the end of each simulation 
cycle. 



D. Lattice geometry, many-body Hamiltonians and 
universal simulations 



In practice, the physical configuration of a device de- 
signed to simulate quantum dynamics is likely to be fixed. 
This may imply that the lattice pattern of the simulating 
qubits can not be changed to adjust the desired simula- 
tion. However, different lattice patterns can be effec- 
tively achieved from a fixed one. For instance, in a two 
dimensional pattern, a system with nearest neighbor in- 
teraction in a triangular configuration can be obtained 
from a rectangular array configuration. This is achieved 
making the subsystems in the rectangular array interact 
not only with their nearest neighbor but also with two of 
their next-to-nearest neighbors in the same diagonal (see 
figure |j). 

As explained in the preceding sections, proposals 
UQSI and UQS2 are best suited to simulating 7V-qubit 
Hamiltonians with one-qubit and two-qubit terms. In 
these cases, the simulation time is proportional to the 
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Figure 4: Illustration how triangular or hexagonal configu- 
rations of atoms with nearest neighbor interactions may be 
simulated in a rectangular array using only nearest neigh- 
bor interactions. Solid lines symbolize interactions between 
atoms. 

simulated time and the control complexity \ is remark- 
ably low. For n > 2, n-body Hamiltonian terms can also 
be simulated in these schemes by exploiting the following 
identity 

e -iH x e e -iH 2 e e iH x o e iH 2 e _ e [H u H 2 }e 2 +0(e 3 ) 

where [A, B] — AB — BA is the commutator of opera- 
tors A and B. A short gate according to a three body 
Hamiltonian H that can be expressed as the commuta- 
tor of two two-body Hamiltonians Hi and H2, that is 
H = —i[Hi,H2], is achieved by concatenating the four 
short two-body gates on the lhs of Eq. (^). As an exam- 
ple, with Hi = 1 ® <j z ® a z and H2 = <t x ®<j x ®1, one can 
obtain an evolution given by Hamiltonian a x ® a y <8> a z . 
A detailed analysis shows, however, that the time T re- 
quired to simulate a three-body interaction by two-body 
interactions grows quadratically in the simulated time T 1 . 
Arbitrary n-body interactions can be simulated in this 
way, with similar polynomial scaling of the simulation 
time. 

Also (i-level systems, d > 2, can be simulated by a 
qubit-based scheme for quantum dynamics simulation. 
For this purpose one needs to group the qubits into sub- 
sets of I neighboring qubits, where I is the smallest integer 
greater than or equal to log 2 d. Single-particle Hamilto- 
nian terms for each d-level systems will correspond to I— 
qubit dynamics. Interactions between n (i-level systems 
will require the simulation of (n x ^)-qubit interactions. 

In this way, a set of qubits with a switchable two-qubit 
interaction and the ability to enforce independent local 
unitaries on each qubit is sufficient to perform univer- 
sal simulation of quantum dynamics, as is the case of our 
second proposal. Instead, when all the control evolutions 
are homogeneous, as in scheme UQS1, only Hamiltonians 
that are invariant under translations can be simulated. In 



the next section we will indicate possible ways to break 
this translational symmetry, thereby endowing scheme 
UQS1 with the capability of performing universal simu- 
lation of quantum dynamics. 



VI. DISCUSSION 

As argued in section I, a general purpose quantum 
computer may only be feasible in several decades from 
now. Instead, a quantum simulator — a device designed 
with the specific purpose of simulating other multi- 
particle quantum systems — could render the simulation 
of quantum dynamics feasible in a much closer future, as 
present experimental results encourage us to believe. 
This paper intended to provide a connection between pre- 
vious theoretical considerations [Q || and the actual re- 
alization of a device for simulating quantum dynamics. 
With this purpose we have put forward two specific ex- 
perimental set-ups and practical mechanisms to attain a 
large class of Hamiltonian evolutions. 

Certainly, much future work is required to complement 
the brief analysis presented here. In this last section we 
would like to describe some of the aspects that should 
be further investigated, and initiate their analysis. In 
particular, we shall examine some of the limitations of 
scheme UQS1 and put forward possible solutions. The 
effect of several kinds of errors is also investigated numer- 
ically in the context of a specific application of a quantum 
simulator of quantum dynamics: the preparation of the 
unknown ground state of a simulated Hamiltonian Hn- 



A. Limitations and imperfections 

Addressability. One of the most significant limitations 
of proposal UQS1 is the lack of independent addressabil- 
ity of individual qubits. In present day experiments per- 
formed with optical lattices, the distance between trap- 
ping sites, d, is smaller than the best achievable focus- 
ing width of lasers beams, w. One typically encounters 
d/w ~ 10 _1 — 10~ 2 , so that several tens of neighbor- 
ing qubits are affected by the same manipulation. Al- 
though homogeneous Hamiltonians, that is, those invari- 
ant under translations, can be still simulated with this 
scheme, it would be most desirable to find mechanisms 
to achieve independent qubit manipulation. In what fol- 
lows we mention three possibilities: 

(i) A spatially dependent magnetic field can be applied 
to shift the energy levels that define the qubits. In this 
way, the original energy difference between the two levels, 
Eq, becomes spatially dependent, E(r), and this allows to 
address the atoms independently. If we shine the system 
with a laser with frequency uol, it will only be resonant 
with the atoms in the lattice site fx, where u>l = 




(ii) The second possibility is to use the spatial depen- 
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dence of the intensity of the lasers. Consider a laser beam 
centered at r (the positions refer to the plane perpen- 
dicular to the laser beam) with an intensity shape pro- 
portional to /(|r*o — ^1), where /(0) = 1. If this laser 
is driving the transition |1) <-* |0) of an atom sitting 
at r, the unitary operation acting on the atom can be 
written as U = exp(— itvof(\?Q — r\)&x)i where vq is the 
coupling constant when the laser is centered at the atom 
position. Suppose we have a chain of atoms at posi- 
tions {rj}, and we want to perform the unitary operation 
V a = cxp(—iTa x ) on atom a leaving the rest unaffected. 
The procedure is to aim one laser to each atom j per- 
forming the unitary transformation Uj — exp(—itjVQa x ) 
for j =/= a and U a = exp(ji a i/o<Tx) on atom a. Since the 
laser beams overlap, the total operation on atom j is 

N 

V s = exp(-iis J2tkf(\rj - r fc |)(-l) 4 *-tr B ). (40) 

k=l 

Now, the condition to be applied to perform the desired 
unitary operation on atom a is Vj — 1 for j ^ a and V a = 
exp(— iTUx) and it is obtained by solving the following 
system of linear equations, 

N 

Ti = ^o5>/(b - r fe |)(-l) 5 -, Tj = r6 ja . (41) 

k=l 

(Hi) The last option is to increase the distance be- 
tween the trapping sites d. This can be accomplished by 
changing the geometry — namely the relative angle of 
the two interfering lasers that form the trap. 

Imperfections. So far in this work we have assumed the 
ideal situation where no errors occur during the manipu- 
lation of the physical set-up that constitute the quantum 
simulator. In practice, however, one has to deal with de- 
viations from this ideal regime, due to several coexisting 
sources of errors: uncontrolled interactions of the atoms 
(ions) with the environment will cause decoherence in the 
A-qubit space; timing errors in laser pulses will trans- 
late into unitary gates (either local or non-local) that 
differ from the intended ones; the adiabaticity condition 
in the performance of short non-local gates in scheme 
UQS1 (UQS2) may not be sufficiently fulfilled, resulting 
in the excitation of the motional degrees of freedom of 
the atoms (ions); 0(t 2 ) corrections in the expansion of 
Eq. ([?]); and others. 

It is beyond the scope and possibilities of the present 
paper to analyze the effect of such imperfections. We 
have nonetheless performed a number of numerical sim- 
ulations involving up to 9 qubits, in the context of a 
particular application of a quantum simulator, namely 
the study of the ground state and first excited states of a 
multi-qubit Hamiltonian. The results, described in the 
next subsection, suggest the scheme is considerably ro- 
bust against several kinds of errors. 

It is also worth mentioning here that decoupling tech- 
niques [p9| can be applied in our schemes to overcome 



some of the difficulties we have mentioned above. In 
particular, undesired local phases appearing during the 
simulation of Hamiltonians can be eliminated using fast 
homogeneous local unitary transformations. This proce- 
dure may also be useful to eliminate inhomogencous lo- 
cal phases originating in (uncontrolled) variations of the 
magnetic field while using method (i) above to enhance 
the addressability in proposal UQS1. The procedure is 
based on the basic observation that a one-qubit evolu- 
tion according to Hamiltonian Hl = o~ z can be undone 
by a similar evolution according to Hamiltonian 

VH L V* = -a z = -H L , (42) 

where {V = iaxjV^} denote fast local unitary opera- 
tions, whereas the two-qubit Hamiltonian Hnl = o~ z ®o~ z 
is left unchanged when the two qubits are similarly ro- 
tated, 

(V®V)H NL (V® V? = H NL . (43) 

Therefore, in proposals UQS1 and UQS2, where the con- 
ditional pushing of atoms or ions produces a short gate 
U according to a Hamiltonian of the form 

E><+E>^\ (44) 

a ab 

a gate according to only two-qubit Hamiltonians H^l 
can be accomplished by performing U, a rotation of each 
qubit according to V ' , another gate U and a final rotation 
of each qubit by V. 



B. Application: studying the ground state of a 
multi-particle Hamiltonian 

The following, final discussion considers a particular 
application of a universal simulator of quantum dynam- 
ics, namely the study of the ground state of a simulated 
Hamiltonian. This is a subject of remarkable interest in, 
for instance, the context of quantum phase transitions 
0. Once the ground state of a Hamiltonian has been 
attained in the A-qubit device, one may attempt to ex- 
perimentally determine its properties, such as correlation 
functions, as well as study the propagation of externally 
induced perturbations. 

The simulation of quantum dynamics by using quan- 
tum systems does not reduce to achieving a controlled 
evolution according to some desired Hamiltonian. Also 
mechanisms to prepare some convenient initial state for 
the ./V-qubits as well as to make measurements that give 
information on the simulated dynamics are needed. In 
section |iy] we already mentioned how to prepare all the 
qubits in the initial state lO) 18 ^, as well as how to per- 
form measurements on the qubits. However, in order to 
study the low temperature dynamics of a given Hamil- 
tonian H j\j , the initial state |0)® Ar must be transformed 



into some other (possibly mixed) state of Hn with much 
lower energy. 

This can be achieved, for instance, by coupling the 
N qubits to a thermal bath at low temperature, fol- 
lowing the line of thought of where conditions on 
the coupling Hamiltonians between the system and the 
bath were studied, so that their joint evolution leads to 
thermal equilibration of the system. In the present case, 
however, with a sequence of gates being continuously ap- 
plied to the A^-qubits, the results of |30| may not apply 
straightforwardly, and a careful study of which kind of 
coupling Hamiltonians lead to thermal equilibrium might 
be required. 

An alternative approach is to make an adiabatic intro- 
duction of the Hamiltonian Hn from some local Hamil- 
tonian whose ground state \^o) can be prepared. 
Suppose that the aim is to bring the N qubits, initially 
in state l^o), into the ground state of H^- This can 
be achieved by simulating a time-dependent Hamiltonian 
H{t) of the form 

H(t) = k(t)H% + (l-k(t))H N , (45) 

where k(t) is a monotonic function that smoothly goes 
from 1 to during a time T s im- The time T s ; m nec- 
essary to obtain the ground state \^>) depends on en- 
ergy difference E e {t) between the ground state and the 
first excited state of H'(t), and must typically fulfill 
Tsim ^ max{ t } E e (i)~ l . Strictly speaking, the ground 
state \^} of Hm will only be achieved in the limit of 
large T S i m . However, for a sufficiently large time T S i m 
compatible with the above condition, the final state of 
the system will be a mixed state p with a large projection 
on the ground state \^>) of Hn and some other contribu- 
tions corresponding to the first excited states. Thus, one 
can view p as a thermal state with very low temperature. 

Numerical simulations of this quantum simulation, for 
up to 9 qubits and including several sources of errors, 
seem to indicate that, indeed, the final mixed state p 
of the A^-qubit system is essentially the ground state of 
-ff/v, together with small contributions of the first excited 
states (see figures (||)-(^)). That is, the final state of one 
such quantum simulations could be used to explore the 
low temperature properties of the simulated Hamiltonian 
Hn- 
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Figure 5: Obtainment of the ground state of a string of 
N = 7 qubits under the presence of dipole-dipole interac- 
tions by starting in the ground state of the initial Hamilto- 
nian H% = ^2i a ai"' Cg) ai a+1 ^ and adiabatically changing the 
simulated Hamiltonian to the dipole-dipole Hamiltonian Ho 
(see section ^) . The simulation includes timing errors in the 
local and two-qubit gates, and also errors due to the first or- 
der approximation in the expansion of Eq. (||) . In Fig. ^(a) 
we assume timing errors in the local laser pulses of 1 percent 
and in the interaction times of 0.5 percent. 100 time steps are 
performed in the adiabatic introduction of Ho, with a value 
of 8i in Eq. ( poj ) of 0.1 at each step. The curve corresponds to 
the fidelity |(*(t)|>I> slm (f))| 2 of the state |* slm (i)} obtained 
through simulation with respect to the instantaneous ground 
state |*(t)) of H(t). Fig. [|b) shows the fidelity of the final 
state with respect to the ground state of Ho for errors of 0, 
1, 2, 3 and 4 percent in both, interaction time and timing of 
local laser pulses. In this case we have 9i = 0.025, and the 
adiabatic change of the Hamiltonian is performed by using 
100 (o), 250 (x), 500 (+) or 1500 (*) time steps. 
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Figure 6: Obtainment of the ground state of Dipole-Dipole 
interaction Hamiltonian Hd (Eq. (j2^)) of a string of TV = 9 
atoms in a ID optical lattice with error parameters of one 
percent for two-qubit gates and timing of laser pulses, and 
6i = 0.025 in Eq. (E0J). The adiabatic change of the Hamil- 
tonian ^ a ai <8> Oi^ 1 ' into Hd is performed by means of 
50 (a), 100 (b), 500 (c) time steps. The histogram shows the 
weight of the state obtained via simulation in the (eventu- 
ally degenerated) eigenspaces Pj corresponding to the sorted 
eigenvalues of Hd, the lowest of which are plotted in Fig. 
6(d). One observes that only the low-energy eigenstates are 
populated, and that the probability of being in the ground 
state converges toward unity if the adiabaticity requirement 
(corresponding to the change of H(t) in Eq. (fta)) is fulfilled, 
that is, if more time steps are performed. 
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